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SUMMARY 


The conservation equations for mass, mean momentum, and turbulent kinetic 
energy for the incompressible turbulent boundary layer have been solved by a 
finite difference procedure. Mathematical models developed by Glushko (1965) 
for the production, dissipation, and diffusion of the turbulent kinetic energy 
in the flat-plate boundary layer have been modified and used to calculate a non- 
equilibrium boundary layer subjected initially to a large adverse pressure gra- 
dient which is followed by a run of constant pressure. 

Comparisons of both mean and fluctuating flow properties have indicated 
generally good agreement between the calculated results and experimental 
measurements of Goldberg’s (1966) . The best overall agreement with data was 
obtained by reducing the scale of turbulence in the outer part of the boundary 
layer to about 70 percent of the flat-plate values as used by Glushko. The cal- 
culations have indicated that further simple modifications to the turbulence 
scale function and to some of the mathematical models for the t;irbulence terms 
should improve the accuracy of predictions for the Goldberg data. It is 
expected that predictions of equal accuracy should be possible for other arbi- 
trary pressure distributions and wall boundary conditions. 

INTRODUCTION 


The basic difficulty in the calculation or analysis of all turbulent flows 
is the problem of how to relate the random fluctuating characteristics of the 
flow to the mean flow properties. In the Reynolds equation.' of mean motion, for 
example, this problem has generally been met by introducing more or less arbi- 
trary assumptions for the Reynolds stress terms (eqs. (l4)« ,17) , Reynolds (1968)). 
In 19 ^ 5 > Prandtl, Nevzgljadov, and Chou reported on independent investigations 
intended to develop a more rigorous approach to this problem. These investiga- 
tions were based on the idea of using independent differential equations to 
describe the dynamics of the correlations for the turbulent velocity fluctua- 
tions. These equations are derived (Hinze (1959) PP- 250-£ 50, for example) from 
the Navier-Stokes equations of motion and contain terms for double velocity 
correlations (Reynolds stress and turbulence kinetic energy terms), triple 
velocity correlations, and correlations of velocity and pressure fluctuations. 

Rotta ( 1951 ) extended the work of Prandtl ( 19 ^ 5 ) , and, in particular, Chou 
(19^5) , in considerable detail based on advances in theory and new data not 
available during the earlier investigations. These methods were further devel- 
oped and applied to various types of simple flows by Emmons ( 1954 ), Townsend 
(1961), Levin (lb 64 ), and Spalding (1967). An integral ’orm of the turbulence 
energy equation was utilized by McDonald (1966) to compute the mean profiles by 
an integral method. ' 

Kovaszany (1967) and Nee (1967) assumed that the efi ective total viscosity 
obeys a "rat** juation" expressing the "conservation" of the total viscosity. 
Harlow (1967) a more formal derivation (based on the equation r ,r the 
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turbulence kinetic energy) of a rate equation for the eddy viscosity and also 
constructed a "transport" equation for the scale of turbulence by analogy with 
Brownian motion. The same authors (Harlow (1968)) have since derived a new 
transport equation for the dissipation function. Since this function depends on 
the scale of turbulence, the new equation together with their previous rate 
equation for the eddy viscosity, were proposed as the basic equations for a 
general method of computing turbulent shear flows. 

In all of the investigations mentioned above, except the integral method of 
McDonald (19 66), it was necessary to make assumptions regarding the relative 
magnitudes of the various fluctuating quantities in order to simplify and obtain 
solutions to the nonlinear equations involved. Consequently, these methods were 
applied mostly to simple flows or portions of simple flows where some terms in 
the correlation equations could be neglected. Hence, the mathematical formula- 
tions of remaining terms representing the fluctuating quantities could not be 
tested for their degree of generality. In the integral method of McDonald 
(1966), it is still not possible to test the detailed spatial variations of 
these formulations and the turbulent diffusion terms drop out completely upon 
integration across a shear layer. 

The answer to these difficulties is, of course, to obtain numerical solu- 
tions of the complete equations with automatic computing machines. Such solu- 
tions have been given recently by Glushko (1965), Bradshaw (1967), and Nash 
(1968). In the methods of Bradshaw (1967) and Nash (1968) the molecular shear 
is neglected and it was therefore necessary to provide the correct wall boundary 
condition by incorporating the "law of the wall" relation between velocity and 
well shear. Glushko, on the other hand, kept all the viscous terns and used 
formulations with the correct limiting form at the wall. It might be expected 
that this latter approach would therefore be somewhat more general than 
Bradshaw's in that sudden changes in wall-boundary conditions could be nego- 
tiated and extension to compressible flows where the "law of the wall" relation 
may not be generally applicable should give better results. The only computa- 
tion presented by Glushko, however, was for the flat plate, where again the 
degree of generality of his assumptions for the fluctuating flow parameters 
could not be determined. 

It is the primary purpose of the present paper, then, to test the method of 
Glushko in a nonequilibrium, adverse pressure gradient flow and to determine 
whether his formulations of the turbulence quantities based primarily on flat- 
plate data result in satisfactory predictions for the mean properties of this 
more complex flow. Since the ultimate success of these methods depends on the 
assumptions relating the fluctuating properties, comparisons of computed values 
of turbulent kinetic energy with experimental data will also be made. 

SYMBOLS 


C constant in dissipation function, equation (5) 

E turbulent kinetic energy profile, e/U ^ 

/ 00 * 

f 2 2 2^ 

e instantaneous value of turbulent kinetic energy, (l/2) (u' + v' + w' j 

e mean value of turbip.ent kinetic energy, (l/2) (u*^ + v'^ + w'^) 

F velocity profile, u/tJ^ 

H function of r or ter in transport functions 

mean scale of turbulence 

n exponent in definition of i], equation (7b) 

Re length Reynolds number, U x/v 
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V 
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turbulent Reynolds number, 


fg l 


constant in transport functions 


transformed normal velocity, equation (8) 
constant in transport functions 

transformed variables, equation ( 7 ) 


constant in transport functions 

turbulence scale function, equation (l) and tabulations, p. 6 


Subscripts 


ave mean or average value 

d dissipation 

5 nominal edge cf boundary layer 

e edge of boundary layer 

o initial condition 

T turbulent 

t transition 


THEORY 


Assumptions for Fluctuation Terms 

The basic physics of the present method (essentially the same as that of 
Glushko ( 1965 )) have been well stated by Reynolds ( 1968 ) . Briefly, the equa- 
tions of continuity, mean momentum, and turbulence energy are solved simultane- 
ously by a finite difference procedure. The expressions for the fluctuating 
quantities in the turbulence energy equation (production, diffusion, and dissi- 
pation of turbulence kinetic energy) as developed by Glushko (19^5) were based 
on the general approach of Rotta (1951) wherein the dissipation and diffusion 

_ \ fp 7 

teims are assumed to be functions of e, Z, and r = * . The form of these 

functions depends primarily on physical and dimensional reasoning. The mean 
scale of turbulence Z was evaluated from flat-plate data for two point 
correlation coefficients of the longitudinal velocity fluctuations, and was 
taken as a "universal" function of the form 

\ - HI) = ">.35 (1) 


Glushko's expression for the Reynolds s cress utilizes a modified form of 
Prandtl's (19^5) eddy viscosity relation and was assumed as, 

Tjp = - P uV = 4 H(r) or (2) 

where 


r x 


Kr) - <£- 



0 < ~ < 0.f5 

O 

0.75 < f- < 1-25 > (3) 

o 

1.25 < — < » 

r 

o J 



The turbulence production term for the problems considered herein is then 


Production = t, 


dU 

T 3y 


= a p 





The dissipation term was written as, 

du. ' du. ' _ _ — 

€ £ v 3x7 S3T- * vC L 1 + H ( CT > a k £I “5 (5! 

J J 

where Z^ = Z, taken as the same function as in equation (l), and H(icr) is the 
same function as given by equation (3) except that r is replaced by icr. 


Glushko assumed that the total diffusion of turbulence energy was due to 
the gradient of e and specified the corresponding diffusion coefficient to be 
the same quantity given in the square brackets of equation ( 5 ) • Hence, he 
obtained for the diffusion terms 



It is obvious that the generality of these assumed expressions for the 
production, dissipation, and diffusion of turbulent kinetic energy can only be 
determined by comparison of final results with data. 


Computing Equations and Procedure 

The computing procedure is essentially the same as that of Glushko ( 1965 ) 
except that the equations are transformed to similarity type coordinates. The 
main reason for transforming to similarity type coordinates is to provide scale 
factors that, in terms of the transformed variables | and tj, reduce or remove 
the rate of increase in boundary- layer thickness with distance | along the 
surface. The number of computing steps. At), required across the boundary layer 
to obtain desired accuracy in the finite difference procedures can thereby be 
reduced and kept more nearly constant. Also, in a region of approximate local 
similarity, the streamwise step size A£ can be increased since, for this situa- 
tion, the rate of change of the dependent variables with § is much reduced. 

The transformed variables are defined as 


r x u . 

i(x) = / ^ dx 

(7a) 

J o 


U oo 


H(x,y) = _ y 

(7b) 


l \ — 7 V \ « / 

v(25) n 

where n is, in general, a variable function of | and is determined from the 
requirement that tj^ (the boundary- layer thickness in the transformed coordi- 
nates) is constant. 


The transformed normal velocity is defined as 

k 2n — « 

x F + ( 2 {) n I- 


V = 


u 


00 


u 


(8) 


and the final computing equations are then written as 
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Momentum.- 


(25) 2 ” P |£ + V || 


/ofc x2n dU 

lil) * (i . p2) 

U &l K } 




«) 


iaulent kinetic energy . - 

+ v 

Sf v 3n 


(25) 2n F ST + v |f =, - ^ FE + (M - 1) (10 


T71\ 2 


+ Si ( D i) ' 


(2|) 2n DE 

„ 2 2 
R 5 q> 


Continuity ♦- 


where 


(a*) 2 ” |f + + (a?) 2 " f (j 


+ In 2|] = 0 


p _ U_ B - e 

F " u «o ' E " it 2 

00 y 


(9) 


( 10 ) 


(11) 


( 12 ) 


M = 1 + or H(r) = 1 + a H(r) cp /E R r 

o 

D = 1 + a K r H(tcr) = 1 + aic H( ter) cp Je R & (13) 

U 5 

R = °o 


The boundary conditions to be applied to the above system of equation are: 

T) = 0: F = E = 0, V = 0, or V = (2£) n X- 

°° (i4) 

tj -»»: F -»1.0, E -» 0, or E -»E K ' 

where E e is one-half the square of the free- stream turbulence intensity. 

The system of equations (9) - (ll) along with auxiliary functions for M, 
D, and n(x), is solved by a linear implicit finite difference procedure. This 
procedure combines certain aspects of the methods given in Glushko (1965) 
Blottner (1964) . The external velocity A and its derivative dU dx must be 
specified functions of x and profiles of V, F, and E must be specified at 
the initial station £ 0 . For details of the procedure and results, the reader 
is referred to a forthcoming NASA publication by the present authors. 


RESULT AND DISCUSSION 


The effect of some modifications to the method of Glushko ( 1965 ) on both 
mean and fluctuating flow properties will be presented for flat-plate flow and 
for one of the experimental flows of Goldberg ( 1966 ) with a large adverse 
pressure gradient. The principal modifications considered are to the l/b 
function (eq. (l)) and to the dissipation and diffusion terms (eqs. (5) and 
(6)). Also, as many of the Stanford cases as possible will be run with the 
"standard" inputs of velocity distribution and initial inputs of Cf, H, 0, and 
velocity profiles. These results will be presented on the standard output plots 
of H, Cp, and Rg. 
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ft should "be noted that in all calculations "by the present method the skin 
friction has been computed from the correct limiting form evaluated at the vail 


as. 



(15) 


The *p functions used in the present calculations vere obtained from the 
follouing tabulated values, vhere the cp „ function is based on the Glushko 
(1965) result. 05 


y £ /8 

’.53 

1/5 

’.25 

^.20 

0 

0 

0 

0 

.2 

.20 

.20 

.20 

.4 

•50 

.25 

.20 

• 5 

• 33 

• 25 

.20 

.6 

•32 

• 25 

.20 

.7 

• 30 

.25 

.20 

.8 

.26 

.20 

.20 

> 1.4 

.01 

.01 

.01 


Linear interpolation between the tabulated values was used in the solutions. 
Prandtl's mixing length relation is thereby recovered in the "law of the wall" 
region (say x g jb < 0.2) where production and turbulent dissipation are approxi- 
mately equal. That is, by equating production and dissipation there is obtained 
(with r » 1) 


1 s -a. I 2 (|nf 

Then with l ~ y, a = 0.2, k = 0.4, C = 3*95 (Glushko (1965)) 


(16) 


r- s (oAy) 2 


5 f 

N 


which corresponds to the Prandtl mixing length relation for turbulent boundary 
layers . 


Flat-Plate Plow 

The calculation was started at Re x = 10 4 with input values of F and V 
from exact numerical solutions to the laminar Blasius flow. The input profile 
for the turbulent kinetic energy was taken as (see Glushko (1965)) 


E 




(17) 


where E * and 


0 ^*/^e are s P ec ^i e 3 constants. Unless otherwise noted, the 


results shown herein were computed with E q * = 2.5 x 10 


-k 


, Tl* . 

and 4 — = 0.4. 


Additional required inputs werej 

U e = 100 ft/sec , v = I.58 x 10" 4 ft 2 / sec, 

.-k 


e ,° 


= 4.95 , 


1 - F = 1 x 10 

e 


■ Si • >• 


4 X 10' 


-if 
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At) = 0.05, 1 - F R = 0.01 and E = 0.82 x 10”®. 

o 5 e,o 

The cp function was used in the flat-plate solution. The values of the 
constants in the transport functions M and D as used by Glushko were 
a = 0.2, re = 0.4, C = 5*93, i* 0 = 110. These values and the function were 

adjusted by Glushko to give agreement with flat-plate flow. In the present 
report, this^ 9.33 function as well as the cp .25 and 9.20 functions as tabu- 
lated on p. 6 , will also be applied to adverse pressure gradient flows. 


Mean flow quantities .- The variation of the form factor H with Re x is 
shown in figure 1. The value of H is at first approximately constant at the 
initial value of 2.592, corresponding to the Blasius solution for laminar flow, 
and then H abruptly decreases at some value of Re x . Increasing the diffus'on 
term by a factor of 3 over the original Glushko form (that is, multiplying 
equation (6) by 3 ) increased the value of x (or Re x ) whore uhe mean profiles 
first began to change from the laminar input shape, by . factor of about I- 1 / 3 . 
(This value of the Reynolds number will be designated ue x ^t, and due to the 
behavior of the mean flow properties can be considered analogous to a transition 
Reynolds number.) However, the H curve appeared to approach the same asymp- 
totic value of approximately 1.4 which is in agreement with the data of 
Wieghardt (1951) considered typical of flat-plate flows. 


The effects of this modification (to the diffusion term) on Cf as 
shown in figure 2 are of the same nature as the effects on H. That is, 
the "transition" Reynolds number is increased by the same factor when the larger 
diffusion term is used but the final asymptotic variation of Cf is in agree- 
ment with that for fully turbulent flow as obtained from Schlicting (i 960 ), 
p. 540, and the data of Wieghardt (1951) • It is also of interest to note that 


V 


was reduced to 1 x 10”® from 2.5 x 10 


when the input disturbance level 
the "transition" Reynolds number was further increased to about 8 x 10 L 


v-4 


The computed values of the mean velocity are plotted in conventional 
profile form in figure 3 and it is seen that the profiles develop from the 
laminar input profile at Re x = 1 x 10^ to turbulent type profiles for 

Re x > 1 x 10 . At this Reynolds number the profile shapes have apparently 
"settled out" to the shape characteristic of turbulent boundary layers as 
indicated by the data of Wieghardt (1951) • When the turbulent diffusion term is 
increased by a factor of 3, the agreement with data is improved, particularly in 
the outer part of the boundary layer. 

The "standard" Wieghardt flat -plate flow (IDENT 1400) was computed with the 
Glushko diffusion (eq. (6)) increased by three, two values of Eq*, and two l/§ 
functions. The 9.25 function appears to give somewhat better agreement with 
data than the 9.33 function, except for the values of Cf at x > 3*0 meters. 


For all the calculations presented herein, the value of n was retained as 
a constant. It was found that appreciable savings in computer time could be 
obtained by the use of a value of n appropriate to the type of flow being 
computed. For example, for the flat-plate problem discussed in this section, if 
the value of n at I = 5 x 10^ was changed from 0.5 to 0.8 the number of At] 
steps required at i = 10° was reduced from 290 to 125* Actually the problem 
had to be reinitialized at £ = 5 X 10 ^ to accommodate the new value of n, but 
the computer time from that point on was reduced by approximately one -half . 


Fluctuating flow quantities . - In figure 4 the computed values of V E 
(total turbulent intensity divided by \j~2) are plotted against y/8p = 0.995 
for several values of Re x and with the same modification to the diffusion 
term as noted previously. The changes in \fiT profiles with Reynolds number 
duplicate those of Glushko with some dependence on the peak turbulent intensity 
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Eo* at the initial station (£ 0 = 4 x 10^). That is, when T^>* is decreased, 
the peak values of \fR at subsequent stations in the transition region are 
decreased, figure 4(b), but the profiles for Re x > 10 ° are of the same 
shape and magnitude regardless of the value of Eo*. This computed behavior 
of the sfE profiles is qualitatively in agreement with e-:perimental observa- 
tions in the transition region. The factor of three increase of the turbulent 
diffusion term improved the agreement with the Klebanoff (1955) data and 
decreased \fW in most of the boundary layer except near the outer edge where 
\TF was increased by factors of two or more. These latter effects would be 
expected due to the gradient type model (eq. ( 6 )) used to formulate the dif- 
fusion term. The simple expedient of increasing this term by a factor of three 
produced better agreement with data for both the \J1T and profiles of r/pe 
(not shown herein). Consequently, the remaining discussion will be limited to 
those results with the turbulent diffusion term multiplied by three. It is 
emphasized that these simple modifications appeared to improve the agreement 
between the calculations and experimental data of the fluctuating and mean char- 
acteristics of the ’'fully" turbulent boundary layer and also improved the agree- 
ment with data of tY a "transition" Reynolds number Rex,t* That is, even the 
minimum transition Reynolds number observed in experiments is still somewhat 
larger than the maximum value obtained here of about 3 x 10 • It therefore 
appears possible that criteria for selecting the "best” models for the turbulent 
terms may depend not only on the results for developed turbulent flow but also 
on the location and behavior of the transitional type flow. 


Nonequilibrium Boundary Layers 

The data from the investigation of Goldberg ( 1966 ) chosen as the test case 
for detailed discussion (his pressure distribution number 3 ) should be a par- 
ticularly sev '"•re test of the method because the boundary layer was first driven 
nearly to separation by a large adverse pressure gradient and then allowed to 
relax toward a flat-plate flow by imposing a constant pressure. Also, hot-wire 
measurements of turbulent shear and longitudinal turbulence intensity were 
available . 


The distribution of external velocity and its derivative with respect to x 
as used to obtain the present results are shown in figure Since there was 

some uncertainty in reading the small graphs published by Goldberg ( 1966 ), two 
alternate velocity distributions are shown in figure 5 (a) and the difference^ 
between them are within the reading accuracy of Goldberg's original figures ar ' 
probably within the experimental errors of the original data. The derivative} 
of the two velocity curves, however, are considerably different as shown in 
figure 5 (b) where dU e /dx as used in the cal culat ions, is plotted against x. 
These differences can become important, as will be shown, when the boundary 
layer approaches separation. Table I lists the values of U/U e and E used 
at the input station corresponding to x = 4 inches (presumably the distance 
from the nose of the 10-inch-diameter test cylinder) . Tne ini tied 5o for this 
station was computed by assuming flat-plate flow at Ue = 85 ft/sec with 
sufficient Reynolds number to give the observed skin friction of 0.00350. The 
initial velocity profile was taken directly from the data plot of Goldberg 
( 1966 ) at x = 4 inches and the initial E profile was taken from the measured 
longitudinal intensity at the same station as 

j (18) 

measured 

since an examination of data where all three components had been measured 
(Klebanoff (1955), for example) showed that this relation wau approximately 
correct for 0.1 < y/5 < 0 . 8 . 


E = 


,2 , .2 . ,2 / , 
u' + v' + w' „ u' 


2U 


00 
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Relation between skin friction and tfW profiles .- The calculated varia- 
tion in skin friction is compared with the experimental data in figure 6. The 
vertical arrows represent the spread in the experimental skin friction data 
obtained by four methods as discussed by Goldberg (1966) . Theoretical results 
are presented for both velocity distributions of figure three values of C 
(dissipation constant), and the three i/5 functions. For velocity distribu- 
tion number 1, C = 3*93> and cp the agreement with data is good except in 
the region of the minimum Cf where the present method overpredicts Cf by as 
much as 100 percent. The agreement xs better when the alternate velocity 
dj stribution number 2 was used indicating the sensitivity of the results in this 
region to the imposed velocity distribution. 

In order to determine the relative importance of the dissipation term (eq. 
(5)) for this particular type of flow, additional solutions with C = 3 and 6 
for pressure distribution number 2 were obtained. An increase in C increases 
the dissipation and reduces the skin friction by an almost constant amount over 
the entire test region and gives improved agreement with data near the minimum 
Cf region at the expense of poorer agreement elsewhere. Since the turbulent 
kinetic energy equation can be paraphrased as 

^ = Production + Diffusion - Dissipation (19) 

it is evident that an increase in dissipation should decrease e, and this, in 
turn, should decrease the turbulent shear which from equation (2) it 
approximately 

T ip - a P \T® * ^ ( 2 °) 

for r > r 0 . Since the skin friction depends directly on the magnitude ol Tm 
in the wall region, the noted decrease in Cf appears reasonable, althc it 

is somewhat surprising to find the almost linear (inverse) relation betwte: the 
magnitudes of the dissipation and Cf. 

A physical "explanation" for the improved agreement beiwet' che computed 
and experimental skin friction in the minimum Cf regj on, as caused by larger 
values of C, is to be found in the behavior of the e profiles which are shown 
in figure 7* 'The theoretical values in is figure were computed with the 
diffusion term taken as three times Glushko diffusion (3 X eq. (6)), velocity 
distributions 1 and 2, different values of C, and different i/5 functions. 

The first thing to note is that all these various modifications had only minor 
effects on the magnitude and distribution of \[E. An increase in C does 
reduce /E, as it Siiould according to equations (19) and (20), but the best 
agreement with the data is generally obtained with C = 3* 93* The input profile 
of \fE (fig. 7(a)) is similar to the flat-plate profiles of figure 5» As the 
minimum Cf region is approached (fig. 7(b)) the peak in both the computed and 
measured \/E profiles moves away from the wall and increases in magnitude, and 
there is a corresponuing large increase in the average turbulent kinetic energy 
(e)ave across the entire boundary layer. That this increase in (e) ave can 
be associated with a decrease in a dissipation length scale Z^ can be seen by 
noting that for large r, the dissipation (eq. (3)) is approximately 

e = «a (e) 5/2 £L (21) 

where Z<j is defined here as a microscale of the turbulence. Then if C is 
regarded as a "universal" constant, the larger values of C as used for the 
solutions of figure 6 should be considered as equivalent to corresponding 
decreases in the square of the dissipation scale, Z^ 2 . If Z is considered 
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an integral scale, the relation between l and for isotropic turbulence 

(Hinze (1959) > P* l85) may be writt * as 



and this relation is seen t® be in agreement with the results of figure 6 as 
related to the change with x in /E profiles of figure 7> since as (e) ave 
increases, l ^ should decrease which, according to the above reasoning, accounts 
for an increase in dissipation, a reduction in /W, and finally, the reduction 
in Cf, as computed . It can then be expected that a functional relation for Z<j 
of the type given by equation (22) would result in good agreement over the 
entire test length for this case. The above explanation depends on the 
assumption that the integral scale Z is not affected by (e) ave . 

Boundary-layer thickness parameters and mean velocity profiles .- The 
computed values of 9 are compared in figure 8 with the experimental data of 
Goldberg (1966) . The differences in the values of 9 between the theoretical 
calculations and the experimental data are small, but, as is well known, these 
small differences can lead to large effects on C^. In general, the velocity 
distribution number 1 gives better agreement with" the data but it is obvious 
that a comparison of P alone is not sufficient to Judge the accuracy of a 
method when C-p -> 0 . 

The shape parameter K is a more sensitive indicator of the accuracy of a 
method as shown by figure 9 where the effects of the two velocity distributions, 
the values of C, and the l/b functions are shown. The velocity distribution 
number 2 and the largest value of C give the best agreement with the experi- 
mental data for z/5 = cp.jj* However, the use of the other two z/5 functions 
gives the best overall agreement with the data, and bracket the data in the 
vicinity of the peak H. The reason for the better agreement of H with data 
for l/b = cp.20 an< * ’ ) .25 is apparent from comparisons of computed velocity 
profiles (not shown herein) with the data. In the vicinity of the minimum Cf 
the agreement between the computed and experimental velocity profiles was poor 
for l/b = 9 ^ regardless of the velocity distribution or the value of C. 
However, wheri cp ^o and vere used, the theoretical results bracketed 

the data. 

Discussion of turbulence scale functions .- It has already been noted from 
figure 7 that ^ is relatively insensitive to changes in free- stream velocity 
distribution, C, and the l/b function. The computed mean velocity profiles 
were insensitive to velocity distribution and C, but could be affected con- 
siderably by a change in l/b. Since the turbulent shear, production, and 
dissipation are the dominant terms in the equations, it is apparent from the 
above discussion and from the form of these terms (see eqs. (2) - (5)) that the 
only way to effect the mean velocity profiles, and hence the H values to any 
appreciable extent is to modify the turbulence scale function l/b. In the 
previous section, a change in C was related to a change in a dissipation 
microscale 2^, but from equation (l6), a change in C corresponds also to a 
change in Prandtl's mixing length constant. Since Prandtl's mixing length rela- 
tion is known to apply even in an r ' -erse pressure gradient in the law of the 
wall region, the l/b function wac not changed for y/b < 0.2. The Goldberg 
(1965) experimental values of mixing lengib ; „ tr ’.hat l/b should be 

decreased in the region cf y/b > 0. T’ decrease would reduce the 

turbulent shear and thereby result i, b rnent with the experimental 

velocity profiles, the l/b function- 1 ged in the manner shown by the 

tabulations on p. 6. The eddy viscosity function e/p does not depend on free- 
stream velocity distribution or C (since \jE* was independent of these parame- 
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ters) but is directly dependent on l/b and agrees better with the data for 
l/b = <p Ibis direct dependence of e/p on l/b is, of course, the reason 

for the’marked effect of the l/b function on both Cj* (fig. c) and u/b'^. 

The success of these modifications to l/b indicates that further minor adjust- 
ments to these functions should give any desired degree of agreement with 
experimental data. In particular, the trends in \[W (fig. 7) for this problem 
indicate the magnitude of the change in l/b should probably depend on the 
level of (e) ave . Spalding ( 1967 a) has applied to free sheen* flows a differen- 
tial equation for the mean scale of turbulence. This equation was based on 
Rotta's hypothesis (1951 , Part II) and relates the scale of turbulence to 
(e) ave and t,j.. It is possible that the use of a similar relation may improve 
the predictions of this present method. 

Other adverse pressure gradient cases .- Additional adverse pressure gra- 
dient cases have been computed with the "standard" velocity distributions as 
supplied for the Stanford cases and for various values of C. The l/b func- 
tion used for these additional cases was generally the 9.25 function as tabu- 
lated on p. 6. Exceptions are noted on the figures. Some difficulty was 
experienced in obtaining these solutions since input values of E were not 
available, except for Bradshaw's (IDEHT 2400). When input E values were not 
available, a trial and error reinitializing procedure to obtain input E was 
used. This procedure was satisfactory in some cases such as Ludvig Tillman 
(IDE 3 TT 1200) but in other cases such as Moses (IDENT 38 OO), the final results 
appeared to be quite sensitive to input E profiles as well as input profiles 
of V. 

CONCLUDING REMARKS 

The equations for the incompressible , turbulent boundary layer with 
constant fluid properties have been solved by a numerical procedure in simi- 
larity type coordinates. Comparisons of calculated values for both mean and 
fluctuating flow properties with experimental measurements in nonequilibrium 
boundary layers as well as the flat-plate boundary layer have indicated 
generally good agreement. 

For the flat-plate calculation, the laminar Blasius velocity profile and 
arbitrary small "disturbance type" profiles for the turbulent kinetic energy 
were used as initial conditions at a Reynolds number of 10 . As the calculation 
proceeded, little change in the mean profiles of velocity and turbulent kinetic 
energy were noted until at seme downstream station, depending on the level of 
the input disturbance and modifications to the turbulence terms, rather abrupt 
changes began and the subsequent mean velocity profiles were qualitatively 
similar to those observed experimental ly in the transition region between 
laminar and fully turbulent flow. The Reynolds number at which these changes in 
the mean profiles were first obtained in the calculation can therefore be termed 
a "transition" Reynolds number and its dependence on the level of the input 
turbulent kinetic energy was shown by Glushko; it is shown herein that modifi- 
cations to the models of the turbulent terms also affected this "transition" 
Reynolds number. 

It was found that when the turbulent diffusion term of Glushko's was 
increased by a factor of 3 > the agreement with experimental values of mean 
velocity and the ratio of turbulent shear to turbulent energy was improved in 
the outer portion of the fully turbulent boundary layer for both the flat-plate 
and nonequilibrium flows. It was also found that when the dissipation term was 
changed for the nonequilibrium flo* , the skin friction was reduced by almost 
constant increments that depended directly on reductions in the square of the 
dissipation scale. Analysis of these results indicated that the micro- 
dissipation) scale may be related to the integral scale of turbulence in about 
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the same way as for isotropic turbulence. The best overall agreement with 
measured values of skin friction, form factor, mean velocity profiles, and 
fluctuating properties was obtained by reducing the value of the turbulence 
scale from a peak of 0.33 to 0.20 or 0.25 of the boundary- layer thickness in 
the outer part of the boundary layer. The linear relation of turbulence scale 
with distance from the wall, in accordance with Prandtl's mixing length theory 
and as used by Glushko, was retained in the "law of the wall region." It is 
concluded that simple modifications to the turbulence scale function and to the 
turbulent fluctuation terms as modeled by Glushko result in accurate predic- 
tions of mean and fluctuating characteristics of turbulent and transitional 
boundary layers with arbitrary boundary conditions. 

TABLE I.- INITIAL VELOCITY AND ENERGY PROFILES AND OTHER 
INPUT DATA FOR ADVERSE PRESSURE GRADIENT CASE 


Values based on data of Goldberg (1966) 



* 8 x 10^ , v = 

1.8 x 10~ k ft 2 / 

1 

sec , p = Oj 

“e, 

= 77*8 ft/sec 

, A$ = 625, ATI 

= 0.025 , n 

1 - 

F = 8 x 10“ 5 
e 

■(ID e = 1 x io ’ 

.5 

’ E e - 1 * 


n 

\ 

\ 




■■■■ 


0 

0 

0 


.1 

.215 

.22 


.2 

.415 

.54 



.5^0 

1.05 


.4 

.580 

1.04 


.5 

.608 

.92 


.6 

.626 

.86 


.8 

.644 

.78 


1.0 

.66 2 

.72 


1.2 

.681 

.675 


1.4 

.700 

.64 


1.8 

.731 

.615 


2.2 

.760 

. 579 


2.6 

•790 

•530 


3-0 

.818 

.474 


4.0 

.880 

•321 


5.0 

.924 

.203 


6.0 

.956 

.108 


6.5 

.968 

.069 


7.0 

.975 

.039 


7.5 

.982 

.010 


8.0 

.991 

.01 


8.5 

.996 

.01 


9.0 

•999 

.01 


.-4 


x 10 


-2 
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Blasius solution 







Figure 2.- Effects of modification to diffusion term and of different values of E q * (maximum 

input E at O on C r distribution for flat-plate flow. 
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Figure 3*- IX /elopment of velocity profiles from laminar input at 
Re = 1 x lCr to turbulent at Re > 1 x 10° for dp/dx = 0. 
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Same as Glushko 

Three times Glushko diffusion 

Data of Klebanoff (1955) 
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Figure 4.- Concluded 




(a) Velocity variation with x. 

Free- stream flow distribution from the Goldberg (1966) pressure distribution 
number 3 as used in the present calculations. 


Velocity distribution 




Velocity 

distribution 



Figure 6.- Comparison of calculated skin friw-ion with measured values of Goldberg (1966) for 
different velocity distributions, values of C, and l/b functions. Diffusion = 3 times 
Glushko diffusion and fi(r) « 1 in outer pare of boundary layer. 





(b) Variation of velocity derivative. 
Figure 5.- Concluded. 




12 inches and input profile used to start solution at x = k inches. 
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Figure 7*- Continued 
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Figure T»- Concluded. 
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